Critical properties and Bose Einstein Condensation in dimer spin systems 
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We analyze the spin relaxation time 1/Ti for a system made of weakly coupled one dimensional 
ladders. This system allows to probe the dimensional crossover between a Luttinger liquid and a 
Bose- Einstein condensate of magnons. We obtain the temperature dependence of 1/Ti in the various 
dimensional regimes, and discuss the experimental consequences. 



Bose-Einstein condensation (BEC) is a particularly 
spectacular manifestation of quantum physics and the in- 
discernibility of the particles. Traditional systems where 
to observe such a phenomenon are bosonic systems such 
as 1 - Helium 4 or more recently cold atomic gases 2 -. How- 
ever BEC is also present in many other condensed mat- 
ter systems. Superconductivity can be viewed as a BEC 
of Cooper pairs and BEC has also long been sought in 
excitonic systems 3 -. Because spin- 1/2 operators can be 
faithfully represented by hard core bosons^, spin systems 
are also very natural systems where to look for such a 
phenomenon 5 -! 6 -! 7 - and various applications have thus been 
investigated. However, in practice, both the experimen- 
tally large exchange constants and, on a more theoretical 
level the hard core constraint, which corresponds to very 
strong repulsion of the bosons, restrict considerably the 
usefulness of these concepts for simple regular spin sys- 
tems. In the extreme case of one dimension it is in fact 
much more fruitful to think of spins-1/2 in terms of spin- 
less fermions^ than in terms of bosons. 

Nevertheless, it has been realized 7 - that the BEC 
physics is extremely useful to describe the spin sys- 
tems made of coupled spin- 1/2 dimers. In such sys- 
tems the application of a magnetic field induces a zero 
temperature quantum phase transition (QPT) between 
a spin-gap phase with zero magnetization and a phase 
where some of the dimers are polarized in a triplet state 
giving a nonzero magnetization. This transition has 
been shown! to be in the universality class of BEC, 
where the bosons represent the dimers in the triplet 
state. The BEC order parameter is then the staggered 
magnetization in the plane orthogonal to the applied 
field. These predictions have been spectacularly con- 
firmed by magnetization 9 - and neutron 1 ^ measurements 
on TlCuCi3 and related dimer compounds. The ac- 
tive discovery and analysis of new compounds in which 
this phenomenon could be observed thus triggered fur- 
ther theoretica l 11 ' 12 ' 13 ' 14 ' 15 ' 16 ' 17 ' 18 ' 19 and experimental 

work 2C i 21 1 22 i 23 1 24 i 25 1 26 1 27^8 1 29^0 Qn ^ magnon BEC . 

In connection with the above mentioned magnon BEC, 
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FIG. 1: Summary of the temperature dependence of 1/TI. 
T c (h) represents the BEC condensation temperature. Ao is 
the gap in the absence of magnetic field. The temperature 
dependence is indicated for the various regimes (see text). 



an important question is to understand how the dimen- 
sional crossover, between the one dimensional and the 
three dimensional behavior can occur in systems made 
of weakly coupled ladders. From a theoretical point of 
view, in such systems, at temperatures much larger than 
the inter-ladder coupling, the system can be viewed as a 
collection of one dimensional ladders, with the very pe- 
culiar associated one dimensional behavior— i 31 ' 32 ' 33 , more 
fermionic in nature, and in the universality class of Lut- 
tinger liquids at high field 8 . At lower temperature, in- 
terladder coupling cannot be ignored anymore, and the 
system falls into the universality class of magnon BEC 
condensates. From the experimental point of view, such 
a change of regime between the one dimensional limit and 
the three dimensional one is relevant for compounds such 
as BPCP 3 - that are a very good realization of weakly 
coupled two legs ladders, with strong dimers along the 
ladder rungs. To probe these different regimes and the 
dimensional crossover, it is crucial to understand the 
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spin dynamics as a function of temperature and mag- 
netic field. An important probe in that respect is pro- 
vided by Nuclear Magnetic Resonance (NMR) measure- 
ments. Quantities such as the relaxation rate 1/T\ have 
been computed in the condensed phase and in the one 
dimensional Luttinger liquid regime 7,32 . However, pre- 
dictions for their behavior close to the BEC quantum 
critical point or in the other one-dimensional regimes are 
still lacking, although some results for two dimensional 
systems have recently appeared 3 ^. 

In this paper we thus calculate the NMR relaxation 
rate \jT\ as a function of temperature and magnetic 
field in a system made of weakly three-dimensionally cou- 
pled ladders, in the one-dimensional and in the three- 
dimensional regime. Our findings are summarized in 
Fig. ED 

We consider a system of coupled one dimensional two 
leg ladders described by the Hamiltonian 

H = J ^ ' S a ,i,l ' S a ,i+e x ,l J-L ^ ' S a ,i,l ' $>a,i,2 
a,i, 1=1,2 a,i 

-h £ S^ + H' (1) 

a,i,l=l,2 

where a is the ladder index, I = 1, 2 denote the two legs 
of the ladder, e x is the direction parallel to the ladder leg 
and ft. the applied magnetic field. H' is the interladder 
exchange term. We take here the simple form for H' : 



important 



7,39,40,41,42 



This allows to use a pseudo spin 



H' — J' S a .i^2 ■ S a , 



i+e y ,l ' 



"J" / , Sg,i,l' Sg,i+e z ,l (2) 



where e y is the direction of the rungs, and e z is the di- 
rection orthogonal to the ladder planes. We take for sim- 
plicity J' = J", generalization is immediate. The NMR 
relaxation rate T\ can be expressed 3 ^ as a function of the 
retarded spin-spin correlation function. We concentrate 
in this paper on the contribution of the transverse local 
spin-spin correlation, 



- = T lim 



lm(S + (q,uj)S (q,u)) 



(3) 



as one can show that the longitudinal contribution 
(S Z S Z ) gives only subdominant corrections as far as the 
temperature dependence is concerned^ 7 .. 

When the temperature is larger than the interladder 
coupling J' the system can be essentially described by a 
set of independent ladders. This one-dimensional regime 
is well understoo d 7 ' 31 ! 32 . For fields ft lower than the spin 
gap of the ladder Aq, and at sufficiently low tempera- 
ture, the system remains in the spin gap state, and \jT\ 
has an activated form 3 ^. For fields ft > Ao , and such 
that ft — Ao T, the system is in a Luttinger liquid 
regimeiS 2 ^ and one finds l/T x ~ T 1 ^ 2 ^ 1 for ft > A , 
with K the Luttinger exponent. Close to the critical 
field the ladder system can be mapped to a two state 
system for which only the singlet and lowest triplet are 



1/2 representation. For high temperatures, T ^> \h— Aq\, 
the system is in a universal quantum critical regime 31 
with a dynamical exponent z — 2. There, 1/T± varies as 
a power law with a universal exponent which can be de- 
termined by a scaling argument. Indeed, scaling requires 
a local spin correlation function of the form T _1 / 2 /(w/T) 
in one dimension. The requirement that / be analytical 
in bj for T > with a vanishing imaginary part for u = 
implies T/u>x1mf(u/T) -> Ct., leading to l/T x ~ T^ 1 / 2 . 
We note that since at the commensurate-incommensurate 
transition induced by the magnetic field, one predicts^^ 
that K = 1, the behavior of 1/T\ is continuous at the 
crossover from the Luttinger liquid to the z — 2 quantum 
critical regime. The exponent simply becomes constant 
once the z = 2 quantum critical regime has been reached. 
Such a behavior is shown in Fig. [TJ 

At low temperature the system enters a three dimen- 
sional regime, where the interladder coupling must be 
considered from the start. Two distinct regimes occur 
depending on the magnetic field. If h — Ao » J' ', then 
the system can be viewed as a collection of one dimen- 
sional ladders, with a fully developed Luttinger regime 
for each ladder. The coupling between the ladders then 
induces a three dimensional ordering for the magnetiza- 
tion in the XY plane, in the same universality class than a 
Bose-Einstein condensation 7 -. The one dimensional fluc- 
tuations strongly affect the three dimensional ordering 
temperature to lead to T BEC = J( J'/ J)W- 2K \ This 
regime can essentially be described by an RPA treatment 
of the one dimensional ladders. If ft. — Ao <C J' then the 
coupling between the ladder is too strong for the ladders 
to have a one dimensional regime. One must then view 
the system as an anisotropic, but three dimensional, spin 
system. Note that this regime is universal and depends 
only on the presence of dimers in a three dimensional en- 
vironment and not on whether a ladder structure exists. 
The results we obtain in this paper for the 3D quantum 
critical regime thus also apply to compounds for which 
Jj_ 3> {J, J') without any specific relation between J 
and J'. A useful representation 7 ^ is to introduce a bo- 
son excitation for each rung in the ladder. The presence 
(resp. absence) of the boson indicates that the rung is 
the triplet (resp. singlet) state. In the regime with three- 
dimensional coupling, and above the BEC condensation 
temperature, we can describe the coupled dimer system 
(TT]) using the following effective Hamiltonian 7 -: 
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h)b{b k 



L 3 ^ 



frfe+A-,^, (4) 



ki,k 2 ,q 



where e(k) = A + is the dispersion of the S z = 1 

gapped magnons, b\ creates a magnon excitation with 
momentum k and g is an effective repulsion between the 
magnons. Since the spin operator S + is proportional to 
the magnon creation operator w, the correlation func- 
tion in ^ can be obtained from the Matsubara Green's 
function of the magnon operator. Above the BEC tern- 
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perature, anomalous propagators are vanishing, and this 
Matsubara Green's function is obtained as a function of 



the magnon self-energy £( k, ad- 
equation. We therefore have: 



-iO) from the usual Dyson 



J 



1 TP d 3 k ImS(fe,w) 
— ~ lim — / - — -r n ■ 

Ti ^uj (2tt) 3 (uj + h-e(k) - Re£(fc, w)) 2 + Im 2 £(fc, w) 



(5) 
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Just above the BEC temperature, the magnon density 
is small, making a low density expansion applicable. At 
first order— only Hartree-Fock diagrams are present, how- 
ever they do not contribute to the imaginary part of 
the self-energy. The first non-trivial contribution to the 
imaginary part of the magnon self-energy comes from the 
second order graphs shown on Fig. [2j In order to evalu- 
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ate the expression ([3]), we use the fact that when u> — > 0, 
TImI](fc, lu)/lo has a finite limit. This implies that the de- 
nominator in Eq. ([5]) is dominated by the contribution 
coming from the real part of the Green's function. Since 
the contribution coming from the denominator is strongly 
peaked near k — when h is near Ao, the dominant con- 
tribution in Eq. ([5]) comes from the small momenta. The 
width in momentum of the expression in the numerator 37 
being of order T 2 / 3 , it is less strongly peaked near 
k = than the contribution from the denominator for 
T — > 0. This justifies the replacement in the numerator, 
of TImE(fc, lo)/oj by its value for k = 0,oj = 0, yielding 
1/Ti - T[m 3 /2(A - /i)] 1 / 2 lim w ^ lml](fc = 0,w)/27ro;, 
with A the effective gap containing the Hartree-Fock 
renormalization of Aq. 
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FIG. 2: Skeleton self-energy graphs at second order. The 
dotted line represents the interaction and the straight line the 
magnon propagator. The magnon propagator is computed 
with the first order self energy. Because of the low density 
limit close to the quantum critical point this expansion is 
well behaved. 



Following the same steps as in Ref. 43|, the graphs of 
Fig. [2] are evaluated in the limit k — > as: 



J 



Im£(fc = 0,u) 



m g 



47T 3 



-(1 



-(A+ei-u)/T 



(1 - e -( A+e i-")/ T )(l - e -( A+e i)/ T ) 
I 



In 1 



a -(ei+A-w)/T 



(6) 



In the vicinity of the quantum critical point, A <C T and 
the integral ([6]) yields: 



I 



m^ 2 g 2 



(A - h) 3 ' 2 (7) 

Let us note that by a power-counting argument^, we ex- 
pect the real part of the ra-th order correction to the self- 
energy to behave as g n T 1+n / 2 . This implies that at low 
temperature the Hartree-Fock contribution dominates so 
that A(T)-h - gT 3 / 2 in the 3D quantum critical regime. 
Therefore, in the 3D quantum critical regime: 

l/Tx-T 3 / 4 . (8) 



In the BEC regime, we find instead that A(T) — h cx (T — 
T B Ec(h)) and \jT x oc (T-T BEC (h)y 3/2 . This behavior 
is valid when T — TBEc(h) ^> |Ao — h\ where Ao is the 
bare gap. This result breaks down when T is sufficiently 
close to TBEc{h). Then, one has the mean-field behavior 
l/T t - {T-TBEcih))- 1 / 2 . Below the BEC condensation 
temperature, due to the presence of a massless Goldstone 
mode associated with phase fluctuations, one can show 7 
that 1/Ti ~ T. At the other extreme, in the gapped 
regime, |A — h\ ^> T and evaluating © yields 1/Ti oc 
7 i 3/2 e -3A/T Q ur fi nc jings are summarized on Fig. Q] fer- 
al! the different regimes. 
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The predictions of Fig. Q] are directly testable on cou- 
pled ladder systems. For systems made of weakly coupled 
ladders such as BPCE 34 , given the accessible magnetic 
fields and with the condition J' -C J, the one dimensional 
regime is probably the easiest to probe. Thus the clear- 
est test is the T -1 / 2 divergence in the quantum critical 
regime. In principle, one could also observe the gradual 
change of the exponent as one moves inside the Luttinger 
regime towards quantum criticality. However, since the 
the exponent varies between 1/2 and 1/3, this change is 
rather small and this effect is likely to be hard to observe. 
Unless one is very close to the quantum critical point the 
pure one dimensional divergence will merge with the di- 
vergence as 1/(T — Tc) 1 / 2 at the finite temperature BEC 
transition due to the three dimensional coupling. For 
systems for which one is able to probe both the one di- 
mensional and the three dimensional regimes, a clear sig- 
nature of the dimensional crossover will be the divergence 
of 1/Ti in the one dimensional regime followed by a drop 
of 1/Ti in the quantum critical 3D regime. Note that 
going towards the BEC regime will reveal a more com- 



plex behavior where the initial drop when entering the 
3D quantum critical regime will be followed by a diver- 
gence at T = Tbec and then by another drop inside the 
BEC phase. Finally more isotropic systems would show 
a broad 3D regime, allowing to test the quantum critical 
prediction T 3/4 in the quantum critical regime or T in 
the BEC one. Unfortunately in TlCuCLj the 3D quan- 
tum critical regime is affected by a simultaneous lattice 
distortion 4 ^ and one must find other candidates to probe 
for this regime. An interesting one might be provided by 
hydrated copper nitrat o 46 ^ 47 . 
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